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Abstract 

Leslie’s perturbative treatment of the direct interaction approximation for shear tur- 
bulence ( Modern Developments in the Theory of Turbulence , 1972) is applied to derive a 
time dependent model for the Reynolds stresses. The stresses are decomposed into tensor 
components which satisfy coupled linear relaxation equations; the present theory therefore 
differs from phenomenological Reynolds stress closures in which the time derivatives of the 
stresses are expressed in terms of the stresses themselves. The theory accounts naturally 
for the time dependence of the Reynolds normal stress ratios in simple shear flow. The 
distortion of wavenumber space by the mean shear plays a crucial role in this theory. 
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I. Introduction 

Time dependent theories of turbulence can be derived by closing correlations in the 
exact stress rate equations in terms of the stresses themselves. 1,2 The present work is based 
instead on Leslie’s observation 3 that the direct interaction approximation (DIA) equations 
for shear flow 4 can be solved by a perturbation series in the strain rate; the goal is to 
apply Leslie’s ideas to single point time dependent turbulence modeling. This application 
of DIA is natural since DIA is a time dependent theory, but since the fully transient time 
evolution of turbulence cannot be described at the single point level, approximations are 
required. Following Leslie and Yoshizawa, 5 we assume that the imposed strain is small 
enough and varies slowly enough to permit treating it as a perturbation of a background 
state of isotropic turbulence. In this limit, single point relaxation approximations 6 can be 
derived by requiring that the model agree with the universal long and short time limits of 
the two point theory. 7 

The model which results has the “viscoelastic” character emphasized by Crow: 8 in 
simple shear flow in which dUi/dxj = SSnSj 2 , the shear stress r is given at short times 
by the rapid distortion theory (RDT) limit 


r = ye k fj s t 1 ) 

where K is the turbulence kinetic energy and t is time; at long times r is given instead by 
the eddy viscosity formula 

K 2 

T = C„ — S (2) 

e 

where e is the dissipation rate. The short and long time limits Eqs. (1) and (2) are also 
satisfied by previous time dependent models including the time dependent shear stress 
model of Hanjalic and Launder 9 


t = -Cr^t+^=KS (3) 

A 15 

by the time dependent eddy viscosity formula of Ref. 6, and by an RDT-based model 
attributable to Townsend, Hunt, and Maxey. 10 

The present model extends these short and long time constraints to second order effects 
including the inequality of normal stresses in simple shear flow. In this case, the short time 
limit was computed by Maxey; 10 the long time limit is the nonlinear eddy viscosity formula 
of Pope 11 (see also Ref. 12, 13 and 14). An interesting property of the Reynolds normal 
stresses is the time dependence of their ratios: defining bij = 2/3 8ij — Tij/2K } the short 
time ratios are 10 

fen : &22 ■ bss = 8 : —13 : 5 (4) 

whereas at long times, very nearly 1,2 

fen : &22 : fe33 = 4 : —3 : —1 (5) 

Previous theories do not give entirely satisfactory accounts of the transition between these 
short and long time ratios. The complete Launder- Reece- Rodi (LRR) stress transport 
model 15 which generalizes Eq. (3) predicts that the long and short time normal stress 
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ratios are equal; the model is calibrated to predict the practically more significant long 
time ratios correctly. More complex second moment closures which are nonlinear in the 
Reynolds stresses could certainly be made consistent with both limits, but the nonlinearity 
is itself problematical, especially when rapid distortion limits are relevant. 1,2 The RDT 
model of Ref. 10 automatically agrees with the short time limit when it is formulated for 
initially isotropic turbulence but the corresponding long time limits are not satisfactory; the 
model is therefore recalibrated by assuming initially anisotropic turbulence. The degree of 
anistropy is chosen to match the long time limits; this adjustment slightly alters the short 
time ratios. Despite some quantitatively successful predictions of certain time dependent 
flows using this model, 16 this account of time dependent turbulence does not preclude 
further investigation. 

Unlike the advanced second moment closures described in Ref. 1 and 2, the present 
model is entirely linear; the complexity required to account for the second order long 
and short time limits enters through a decomposition of the Reynolds stress into tensor 
components with distinct time dependent behavior. The decomposition is determined 
by a geometric feature of shear turbulence, the distortion of wavenumber space by the 
mean shear. This effect is often treated as a weak mechanism of energy transfer by eddy 
distortion; 18 it proves to be negligible at short times but sufficiently important at long times 
to bring about the transition from Eq. (4) to Eq. (5). The introduction of quantities with 
different time dependence is generally consistent with Weinstock’s suggestions 19 that the 
shear and normal stresses have distinct relaxational behavior. 

It must be emphasized that the present theory is limited to small strain rates with 
rapid distortion theory describing only the short time limit. The more difficult problem of 
turbulence at high strain rates is left to further investigations. 

II. Simplified DIA analysis of shear turbulence 

The DIA generalized Langevin model 20,21 for isotropic turbulence is 

d [* 

—u z (k,t) + J ds rj(k,t,s)ui(k,s) = (6) 

where rj is a damping function and / is a random force. In Eq. (6), the viscous term has 
been dropped. Define the response function G as the inverse operator for the left side of 
Eq. (6) so that 

^G(k,t,s) + j dr rj(k,t,r) G(k,r,s) = 0 (7) 

and 

G(k 1 t 1 t) = 1, G(k 1 t 1 s) = 0 for t < s (8) 

Define the two time correlation function by 

Qij(k,t,s)8(k + k') = < ui(k,t)uj(k' ,s) > (9) 

so that in isotropic turbulence 


Qij(k,t,s) = Q(k,t,s)Pij(k) 
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where k 2 = k • k and Pij is the transverse projection operator 


kikjk 

Then the damping function r] is related to G and Q by 20,21 


ri(k,t,s) 


k=p+q 


dpdq P tmn (k)P mst (p)P ns (q)G(p, t, s)Q(q , t, s) 


( 10 ) 


where 

Pimn( k) = k m Pi n ( k) T k n Pi m ( k) 

The random force fi in Eq. (6) is defined in terms of independent Gaussian random fields 
such that 


• fdk./iGik'.,, > ■ fOk./if'ik'.,, > = Chpk./.^ik • k'l (11) 


by 


/;(M) = P imn ( k) 


dp dq Cm(p ,t)^' n (q,t) 


k=p+q 


( 12 ) 


The DIA equations 20,21 for G and Q follow from Eqs. (6)-(12). The existence of this 
Langevin model insures the realizability of DIA and also brings forward some analogies to 
statistical physics, where Eq. (6) arises as a generic model. 22 

The DIA theory of shear turbulence in the special case of spatially constant mean 
velocity gradient 3,4 can also be formulated as a generalized Langevin model 


d_ 

di 


ui(k,t) + 


f ds ?^j m (k,t,s) u m (k, s) 
o 


S/m (k, t)u m (k,t) + fi(k,t) 


(13) 


where the mean shear operator 18 is 

S im(t) = + 2 k^ktkp^^-it) + 8 irn k s ^^-{t)-^r 

dx m dx m dx r ok. 

and the damping, described by the tensor 


( 14 ) 


?7 s) 


k=p+q 


dpdc[P ip n (^k^P r srn{p)Gp r (^p^t^ S^Q ns (^C[^t^ 5 ) 


( 15 ) 


need not be isotropic. The first two terms in the shear operator of Eq. (14) are the 
production and rapid pressure strain terms; the third is generally described as a relatively 
weak mechanism of energy transfer due to eddy distortion by the mean shear. 18 This 
effect can also be considered purely geometrically as a distortion of wavevector space. The 
random force is defined as above in Eqs. (11), (12); however, the velocity correlation in 
Eq. (11) must not be assumed isotropic. The shear flow DIA equations 4 for Gij and Qij 
follow from Eqs. (13)-(15). 
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The equations of shear flow DIA seem to defy analysis; even numerical investigation is 
hampered by the need to resolve a three dimensional anisotropic energy spectrum. 23 How- 
ever, a useful simplification was suggested by Leslie: 3 introduce an isotropic background 
field satisfying Eq.(6), which we rewrite as 


^u\°\k,t) + f ds ? 7 (0) (A:,t,s)^^ 0) (k,s) 


Ao) 


to emphasize that the damping and forcing are isotropic. Formulated as a Langevin equa- 
tion, Leslie’s model is 


Jb,(k,tf) + J ds r/ (0) (A-,ps) u,( k,s) = /f^k,*) + S, p (k, t)u p {k, t) 


( 16 ) 


in which the damping and forcing are isotropic and independent of the mean shear. Eq. 
(16) can be rewritten in terms of the isotropic Green’s function G^°\ the solution of Eqs. 
(7), (8) as 


u z (k,t) 


ds G m (k,t,s ) [/, l0l (k, s ) + P.,(k)S„(k, s )i. p (k, s )] 


(U) 


This model ignores any effects of the mean shear on the nonlinear transfer mechanism, 
modeled by the damping and forcing in the Langevin equation. This is certainly an over- 
simplification. However, a large scale shear acting against an isotropic transfer mechanism 
provides a plausible picture of shear turbulence. Moreover, the Kolmogorov time scale 
~ A: -2 / 3 is much smaller than the time scale of the mean shear at sufficiently small scales. 
It is reasonable to expect the perturbative treatment to be valid at such scales. 

Because of this simplification, it is important to compare Leslie’s model with the full 
DIA shear flow equations. Denote the solution of Eq. (16) by . Define the first order 
correction to the Green’s function G^} so that 

L J 


,[ i 1] (M) = / ds 


ip 


(18) 


The field can be used to construct a corrected correlation tensor 

Q^(k, f, s)6(k + k') = < «W(k, f)'uW(k / , s) > 


(19) 


Substituting the corrected quantities G^J in Eq. (15) yields a corrected damping 
rjfj and substituting in Eqs. (11), (12) yields a corrected forcing . This approximation 
scheme can be iterated. Define for n > 2 as the solution of 

d 


dt 


Y ] + ds 1] 4 n] = S ip u [ ™ ] + fl n 1] 
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Define and given by analogy to Eqs. (18), (19). If the approach a limit, 

J ... \ n ] \ n ] 

the corresponding limits of G 1 - 1 and Q 1 - 1 are a solution of shear flow DIA. Leslie’s model is 
the first step in this iterative construction. It is therefore a rational approximation which 
can be systematically corrected by further iteration. 

III. Time dependent linear eddy viscosity 

To derive the time dependent linear eddy viscosity, we follow Leslie 3 and expand Eq. 
(17) in powers of the mean strain rate about an isotropic background state u 

u = u (0) +u (1) + ••• (20) 

The random force / which maintains the background field 3 is taken to be of the same 
order as u therefore, t/ 1 -* is given by 

u\ 1 \k,t) = [ ^G^,M)P im (k) S m „ (k,.) 4°) (k,.) (21) 

Jo 

If, corresponding to Eq. (20) the single time correlation function is expanded as 

Q(k,t,t) = Q (0) (k,t) + Q (1) (k,t) + • • • (22) 

where 

Q- J 1) (k,t)<5'(k + k') = < u*p (k, t) u^\k' } t) + u- 0) (k, ^^(k', t) > 
and the higher order correlations are defined similarly, then Eq. (21) implies 3 

Qy’tM) = J <M + 2M' P v 2 Aq etpcM) + (u) 

(23 ) 

where (ij ) denotes index interchange in the immediately preceding term. The k derivative 
terms represent the distortion of wavevector space by the mean shear. A decomposition of 
the Reynolds stress follows from Eq. (22): 

T *j = ~ < u * u j > = T if + T iP + ■■■ 

where 

A'W = -jdkQ\f(k,t) (24) 

Note that is homogeneous of degree n in the mean velocity gradient. In applying 

Leslie’s theory, it will be convenient to invoke the Markovianized version of DIA described 
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in Appendix II. Evaluation of the angular integrals in Eqs. (23), (24) using Eq. (70) leads 
to 3 


pt pOO a r) 7 

r ( i ]\t)= J ds J dk {— G(k,t,s) 2 - — G(k,t,s)k — G(k,t,s)j E(k,s) Sij(s) (25) 
where E(k,s ) is the energy spectrum at time s, and 

_ dlh dUj_ 

13 dxj dx t 


is the strain rate. In Eq. (25) and what follows, the superscript (0) on the background 
field descriptors has been dropped. 

An equation similar to Eq. (25) has recently been derived by Woodruff 24 by a scale 
separation expansion similar to Yoshizawa’s and has been developed into a subgrid scale 
stress model. Our goal is instead the development of a single point model. Eq. (25) 
suggests a decomposition of the first order stress, 


(!) = i T (°) , 1 T (!) 

ij 15 ± ij ' 15 ± ij 


( 26 ) 


where 



t pOO 

ds / dk G(k,t,s ) 2 E(k,s ) Sij(s) 

0 Jo 

r f r°° d 

ds dk G(k,t, s)k [— G(k,t,s)\ E(k, s)Sij(s) 

1 o Jo dk 


( 27 ) 


The term T^ 1 -* exhibits the distortion of wavevector space by the mean shear. The time 
dependent model provided by Eqs. (26), (27) is purely formal because the time evolution 
of G and E is not known but must be found in general by solving DIA to determine the 
fully transient spectral dynamics. At this level of generality, single point modeling is not 
possible; extracting a practical single point model requires further restrictions on the time 
dependence. We can retain some of the time dependence of Eqs. (26), (27) following Smith 
and Yakhot 7 by requiring that the single point model agree with the essentially universal 
short and long time limits of Eq. (27). It must be stressed that this model is necessarily 
of quite restricted applicability; it cannot describe arbitrary complex strain rate histories. 

At short times, we have 

E(k,t) = E(k) + 0(t) 

Sij(t) = Sij + O(t) 


and in view of Eqs. (7), (8), 


G(k, t, s) = 1 + 0(t — s) 2 
^-(M,s) = 0(t - s) 2 
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consequently, 


(28) 


Ti 0| = A-(S„+O(( 2 ), T</>=0(i 2 ) 

Substituting in Eq. (26), there results 

T , 1 , 11 = + ( 29 ) 

in agreement with Crow’s 8 RDT calculation. The long time limit is a hypothetical steady 
state defined by setting E(k } t) and S(t) to constants, setting G equal to its time stationary 
form 

G(k 1 t 1 s) = G(k,t — s) = exp [—r](k)(t — s)] (30) 

from Appendix II, Eq. (68), and setting t = oo in Eq. (27). Exactly homogeneous shear 
turbulence is never naturally in a steady state; this long time limit occurs only locally 
in inhomogeneous flows like the log layer where it can be maintained by strong turbulent 
diffusion. 

After substituting in Eq. (26), the result is the eddy viscosity representation 

T iP = (31) 


where 


4 7°° 2 7°° 

vt = — / dk 9 0 E + — / dk kr]'8 1 E(k ) 
15 Jo 15 Jo 


and 6i are time moments of the stationary Green’s function of Eq. (30), 


( 32 ) 


9 n = 


da a n G(k } a) 2 = 


n\ 


( 27 ) 


n + 1 


( 33 ) 


Results like Eq. (32) which express turbulent transport coefficients in terms of integrals of 
inertial range quantities also arise in Yoshizawa’s theories. 5 In view of the appearance of 
the time integrals 6 n of response functions defined in Eq. (33), Eq. (32) can be compared 
to Kubo formulas 25 expressing molecular transport coefficients in terms of integrals of 
equilibrium correlation functions. 

It remains to construct a relaxation model for T^ 0 -* and T^ 1 -* which is consistent with 
these limits. In both limits, the mean shear can be considered constant; therefore, it 
suffices to consider the time evolution of the quantities 


f (0) = dk ds G(k,t,s) 2 E(k,s) 

Jo Jo 

/»o© pt 

t (11 = / dk / ds kt/ik) (t - s) G(k,t,sf E(k,s) 


( 34 ) 

( 35 ) 


In Eq. (35), the k derivative of Eq. (27) has been evaluated using the time stationary 
form for G of Eq. (30). Approximate the time evolution of t by the relaxation model 


i (0) = k - c^et^ 


(36) 
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where 




is a frequency of the largest turbulent scales and the relaxation constant is to be 

determined. The relaxation model Eq. (36) has been chosen to agree with the short time 
limit of Eq. (34) 

f (0) = Kt + 0(t 2 ) for t — >■ 0 

To match long time limits, set i = 0 in Eq. (36) and replace t ^ by its long time limit 
computed from Eq. (34). The result is 


C't 0) © = 


r m e 

J°° dk e,,E 


(37) 


This procedure can also be applied to the term in Eq. (35). The relation 




dk 


ds [krj' (k) + 2r](k)kr]' (k)(t — s)}G(k,t, s) 2 E(k, s) (38) 


Jo Jo 

suggests the relaxation approximation 

IS 1} = ci^QtS^ - C^ 1} 0t (1) 


(39) 


This relaxation model has the same short time limit, = 0(t 2 ) as Eq. (35). Since only 
the 0{t) term in Eq. (34) has been matched, it would be inconsistent to match terms of 

order t 2 in Eqs. (35) and (39). Instead, is determined by equating the long time limit 
of C , £ 1 ' ) 0f(°) to the long time limit of the first term on the right side of Eq. (38): 


(i)„ f 0 °° dk W d 0 E 

J™ dk 6 0 E 


CgQ = 


(40) 


Now can be determined by setting = 0 in Eq. (39), replacing by its long time 
limit evaluated from Eq. (35). The result is 


(i) Q _ fo°°dk 2 r] kr] 1 Q 1 E 

J^dkkri'OiE 


Cr& = 


(41) 


Substituting the relaxation approximations Eqs. (36), (39) for into the definitions of 
T (*) 

in Eqs. (26), (27), leads to the relaxation approximation for shear stress 


( 1 ) = ± T (o) , 

iQ ± ij ' \§ ±% i 

ij? = ksh - 

if = CfST ™ - C^0T^ 


( 42 ) 



Eqs. (37), (40), and (41) express the relaxation constants C^\C^\C^ in terms 
of the energy spectrum E(k) and damping function r/(k) which can be assumed to have 
Kolmogorov inertial range forms 

E(k) = C K e 2/3 k - 5/ 3 r](k) = C T £ 1/3 k 2/ 3 for k > k 0 

where Ck and Ct are universal inertial range constants, and ko is a (nonuniversal) cutoff 
for inertial range scaling in any particular problem. The dependence of the relaxation 
constants on the cutoff scale ko and on the form of the spectrum and relaxation function 
in the far infrared region k < ko must be investigated. Rewrite Eq. (37), for example, as 

a formula for using the definition of 0: 

(., l [J™dkE(k )] 2 
R £■ /y dk OoE(k) 

and substitute the simplest expressions for E and rj in which E(k) = r/(k ) = 0 for k < ko- 

Then both numerator and denominator scale as k 0 4,/3 ; is therefore independent of 

the cutoff scale ko- Consequently, the relaxation constants can be evaluated in the infinite 
Reynolds number limit in which ko —■ ► 0 and treated as universal properties of the inertial 
range alone. Thus, 

c i 0) = QCkCt 

This result is approximate, but the exact value given by Eq. (37) depends only weakly on 
ko and the particular form of the spectrum and time scale when k < ko- 

IV. Time dependent nonlinear eddy viscosity 


Calculation of Leslie’s expansion to second order is tedious but straightforward. The 
result is 


n (2) (kt)- V I (N) \a (N) — (s)^-(r) + b (N) — (s)^-(r) + c (N) ^( s)^(r) 
2^ 1 \- a dx p [) dx p {)+ ^- [)+ 


1<N <6 


dx p y dxj 


+ d<*(#(r) + + /<">* As)^(r)] 


dxi J dxj 

+ (* 7 ) 

in which the are integral operators 


dx q J dx p 


dx q J dx q 


( 43 ) 


/<■> = / ds G(k,t,s) / dr G(k, s,r)G(k ,t,r)Q(k,r) 

Jo Jo 

I"*- 2 ) = [ ds G(k,t,s) [ dr G(k,s,r)k — [G(k,t,r)Q(k,r)] 

Jo Jo dk 

E r s d 2 

1^ = / ds G(k,t,s) / dr G(k, s, r)k 2 —— [G(k, f, r)Q(k, r)] 
Jo Jo dk 
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( 44 ) 


I il> = f dsG(k,t,s)f dr G(k,t,r)G(k,\ s,r \)Q(k,r) 

Jo Jo 

= [ ds G(k,t,s) [ dr G(k,t,r)k — [G(k,\ s,r \)Q(k,r)] 

Jo Jo dk 

ft ft J2 

/ (6) = / ds G(k,t,s) / dr G(k } t,r)k 2 — [G(k,\ s,r \)Q(k,r)\ 

Jo Jo dk 

understood to act on products of mean velocity gradients by integration over r and s 
and (ij) denotes index interchange. In Eq. (43), are the following geometric 

constants: 



(1) 

(2) 

( 3 ) 

( 4 ) 

(5) 

(6) 

105a( w ) 

27 

-1 

-2 

19 

2 

1 

2 

-1 

1056^ 

20 

6 

-2 

6 

3 

-1 

105c( n ) 

-15 

-15 

-2 

5 

2 

15 

2 

-1 

105d( w ) 

20 

-8 

-2 

-22 

-4 

-1 

105e( N ) 

10 

24 

6 

3 

12 

3 

105 

-4 

24 

6 

10 

12 

3 


and 

G(k } | s,r |) = G(k 1 s 1 r ) + G(k 1 r 1 s) 

As in all derivations of this type 10,12 ’ 13 , 14 ^ these constants arise from integrating even order 
products kikj, .... over spheres k = constant. It is easily verified that these results reduce 
to Maxey’s 10 second order RDT calculation when G = 1. 

In extending the relaxation approximation to second order terms in the strain, Q and 
dUi/dxj can be taken independent of time in Eq. (43). In this case, the integrals in Eq. 
(44) satisfy 

j( 4 ) = 2 1*' 1 '* I*' 5 ) = 2 1*- 2 ) I ^ = 2 1*- 3 ) 


Taking traces in Eqs. (43) and (44), we find that 


dk Q tt (k) 


± Kt 2 dU^dU^ 

15 dx q dx q 


(45) 


exactly equals the energy produced at short times by the shear given by Eq. (29). This 
suggests absorbing the trace of Eq. (43) in the energy equation and retaining only deviators 
in the stress model. Although this procedure is strictly valid only at short times, we will 
generalize it to arbitrary times, recognizing that this is an additional assumption. Another 
important property of the geometric coefficients in Eq. (43) is 


2 (a (JV) + 2a (JV+3) ) 


dU t dUg 
dx p dx p 


+ ( b ^ + 2 + c ^ + 2c( jv+3 ))( 


dUi dU n dUi dU n 


dx p dxj 


dx p dxi 


+ 2(d( JV ) + 2d( JV + 3 ))?^^ 

dxi ox j 

= (a'"' + 2 a '™>)( + ^s ri) + (d<"> + 2<(' «+=>)( + ^S r .) (46) 


dx. 


dx; 


dx. 
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The importance of this condition, which states that the stress vanishes when the strain 
rate tensor Sij vanishes, was first clarified by Speziale. 12,13,14 It permits a consolidation 
of the perturbation expansion by replacing time integrals of S by stresses; compare in this 
respect the procedures of Ref. 26 and 27. These relations have recently been verified in 
Yoshizawa’s theory. 28 

The second order analogs of the quantities of Eqs. (34), (35) are 




dk 


O© 

dk 


d 

dk 






n 


dk 


f 

o 


r 

0 


ds ( krj'y {t — s) l G(k,t,s) 2 
ds ( krj'y (t - s) l G(k,t,s) 2 


dr (krj'y (s - r) j G(k,s,r) 2 E(k) 
dr (krj'y (s — r) J G(k, s,r) 2 X 


(47) 

f ds {krj'y (t — sY G(k ,t, s) 2 [ dr (krj'y (s — rji G(k, s,r) 2 ——(k 2 E) 
o Jo dk 


where 0 < i } j <2 and 0 < i + j <2. The term s^ 0,0 -* is of order t 2 at short times. The 
terms with (i } j) Y (0,0) are evidently of order t 3 or higher; integration by parts 

shows that s^ 0,0 -* and s^ 0,0 -* are also of order at least t 3 . Thus, short time limit of the 
second order quantities is determined by terms containing s^ 0,0 -*. It is shown in Appendix 
I that that these terms can be written as 




(48) 


where the quantities cr^ 1 ), cr^ 2 \ satisfy relaxation equations 


4 11 = Ty°> - c^e^ 11 


A 1 = T. 

•( 6 ) 


a 




= T 


( 2 , 0 ) 

ij 

( 3 , 0 ) 


- c 


(0,0) 

R 


0(7 




- C4 O,O) 0<7 


( 2 ) 

ij 

( 6 ) 

ij 


(49) 


The relaxation constant Cr’°^ is defined by a quotient analogous to Eq. (37) given in 
Appendix I, and the are defined by 


7 -'(ij) 

- L mn 


= (a«+2a< !+3 )){ 
+ {d {l) +2d ( * +3) ){ 


dlJn 

dx f 

dU p 

dXm 


-T (j) + 

pn 1 

T (j) + 

pn 1 


^n T (i) 

dx n pm 


dll* 


j'(j) 


dx„ pm 


(50) 


in terms of the quantities defined in Eq. (27). Absorbing the r integration this way 
is justified by the property Eq. (46) of the second order geometric constants, and taking 
deviators is justified by the discussion following Eq. (45). 

For simple shear flow, the model defined by Eq. (48) predicts that the short time 
dimensionless Reynolds normal stresses deviators are 




2.2 


3 105 


-S 2 t 


, 2 -13 9 2 

622 — S t 

3 105 


2 5 9 9 

= 3106* V 
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in agreement with Maxey’s RDT calculation; 10 in particular, the normal stress ratios are 
given at short time by Eq. (4). However, this model predicts the same normal stress 
ratios at long and short times. The long time ratios can be corrected by adding terms 
to the model which vanish to higher order at short times. As noted earlier, such terms 
originate with the k derivatives in the shear operator Eq. (14). A complete model based 
on analysis of all the terms is presented in Appendix I. Since it is only remotely possible 
that so elaborate a model could be useful in practice, we will present a simple model with 
the correct long time stress ratios. Using the notation for components of the second order 
stress in Appendix 1, this model is 

= A 11 - 3 A’ + 12 A’ - A ) 1 (si) 


where the stress components satisfy relaxation equations Eq. (49) above and 


A 81 = U J) - cT’qA 1 


( 52 ) 


Adding the term — 4<4 8 ) corrects the long time normal stress ratios. Here, the coefficient 
—4 has been adjusted so that assuming 

y'~l(0,0) ^~y(0,l) sy* 

“4 - 


the long time limit is 

1 16 ,SK,2 


6u = 


C* R C° R 105 v e 


1 — 12 2 1 — 4 ,SK o 

v / 22 q* qo 105 e 33 1 nc; ' 2 


C* R C° R 105 v £ 


Setting the phenomenological constant 


C* R = 2.00 

gives 6n = .19 in equilibrium flows for which SK/e = 3.3. As the normal stresses are also 
in the ratio Eq. (5), the long time limit is in good agreement with experimental data. 1,2 

The introduction of different relaxation constants C^ R and C^ 1 ^)r for the shear and 
normal stresses has been advocated previously by Weinstock. 19 Many of the phenomena 
discussed in Ref. 19 could be described by models of this type. 

V. Conclusions 

Leslie’s perturbative treatment of the DIA equations for shear flow has led to a model 
of second order effects with the following properties: 

1. The model consists of linear relaxation equations. This attribute should be compared 
with models nonlinear in the stresses recently criticized by Speziale: 1,2 such nonlin- 
earity can be problematical in connection with rapid distortion limits. 

2. The model correctly accounts for the change in Reynolds normal stress ratio from 
short to long times. This property is interesting from a theoretical standpoint despite 
the limited practical importance of the short time normal stress ratios. 


12 



3. In this model, the shear stress determines the normal stresses, which occur as effects 
of higher order. This property contrasts with stress transport models like LRR in 
which the normal stresses influence the shear stresses. 

The limitation of the model to small strain rates and to strain rate histories which 
are not extremely complex must be noted. This model can be applied to predict normal 
stresses in oscillating flows. This problem has been successfully treated using the RDT 
based model by Mankbadi and Brereton, 16,17 however, as noted previously their model 
slightly modifies the short time ratios. 


Appendix I. Derivation of second order relaxation model 

The second order calculation requires the quantity 

/»OG ft 

t {2) = dk ds (kf]') 2 (t - s) 2 G(k,t,s) 2 E 

Jo Jo 

analogous to and in Eqs. (34), (35). The relaxation model 


t' (2) = 2C ( s 2) et (1) - C^ 2) 0t (2) 


is derived like the model for t^: the result is 




f^dk (krj'fdiE 
j^dk krj'OiE 


C ( r® 


f™dk 2r] ( kr f) 2 6 2 E 

J™dk {k^) 2 Q 2 E 


(53) 


The definitions Eq. (50) of now apply when j = 2. 

Appropriate relaxation approximations for the time evolution of the Eq. (47) are 

i(0’0) =t(°) _ C (°’0) es (0,0) 

■(1,0) = C <A-O) 0S (O,O) _ C'U.O) 0S (1,O) 

^(0,1) =t (l) _ C (0,l) Qs (0,l) 

i( 2, °) = 24 2 ,O) 0 ^ 1 ’°) - C^ 2 ’ o) 0,s (2,o) 

i (1,1) = ddW 0 ’ 1 ) - C'J’ 1) 0s (1 ’ 1) 

^(0,2) =t (2) _ C (O,2) 0s (o,2) 


where the relaxation constants are determined by 

r (o,o)^ _ /o°° dk 2 r\Q x E 
R / 0 °° dk 6iE 

r (i,o) 0 Jo °° dk kri'QiE (i,o) / 0 °° dk 2r 1 kr 1 '6 2 E 

S J^dkB^ R J^° dk kr]'0 2 E 

r (2,o) n _ Io 0 °dk (kri') 2 e 2 E ( 2 ,o) _ /o°° dk 2r 1 (kr 1 ' f0,E 

s dk krj'6 2 E R / 0 °° dk (kr]') 2 6 3 E 

( 1 , 1 ) f 0 °° dHWfe.E dk 2r 1 (kr 1 'f0,E 

s jf 0 °° dk krj'9 2 E R / 0 °° dk (kr]') 2 6 3 E 
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The remaining relaxation constants satisfy 


^(0,1) _ r (0,l) r (0,2) _ r (l,l) _ r (2,0) 

L- D L- 70 L- 70 L- D L. D 

Jri it it it it 


The second order part of the stress is the sum of contributions from + 2 /( w + 3 ) 
for N = 1,2,3. The contribution from + 21 is 


(i) _ 


cr v ' = 


dk 


dsG(k,t,s) 2 / dr G(£g,s,r) 2 £(£;){(a (1) + 2a (4) )(WS + SW T ) 


+ (d (1) + 2d (4) (W T S + SW)} D 

where matrix notation is used and the subscript D denotes deviatoric part. The relaxation 
approximation is 

<j ( i) = T (1 ’ 0) - C4 0,0) <7 (1) (55) 

The contribution from I ^ + 2 i^ 5 -* can be written 


f s d 

dk / ds G(k,t,s) / dr G(&, s, r)& — [(?(&, A, r)Q(&)] x 

do Jo dA: 

IV 2 ) + a^)(VUS + SW T ) + (d< 2 > + d^)(VU T S + SW)} D 
= -3u( 2 ) - cr(3) _ ^(4) + u(5) 


where 

a (2) = s ( 0 , 0 )|( a ( 2 ) + 2a^)[VUS + SVU T ] + (d< 2 > + 2d^)[VU T S + SW]} D 
<t ( 3) = s (1 ’ 0) {(a (2) + 2a (5) )[V17S' + SVU T ] + (d (2) + 2d (5) )[W T S + SW]} D 
<t ( 4) = ^’^{(a^ + 2a ( ' 5 ' > )[VU S + SW T ] + (d< 2 > + 2d^)[W T S + SVU}} D 
a (5) = s (o,oy^ a (2) + 2a ( 5 ))[V!75 + SVU T ] + (d (2) + 2d (5) )[V17 T S' + SW]} D 


and 

dk G(k , f, s) 2 

In the relaxation approximation, 



G(Ag s, r) 2 


d 

dk 


(kE) 


d (2) = T (2,0) - C^°’ o) 0<j ( 2) 
d (3) = cf ’ o) 0T (2 ’ o) - C'J’ o) 0<7 (3) 

(j(4) = t ( 2,1) - c4 o,1) 0<7 (4) 

d (5) = 2T (2 ’ 4) - C'2 ,,O), 0 <t (5) 

where 

(0;0) , _ J 0 °°dk 2r ] d 1 (kE)' 

R j^dke^kEy 


(56) 


(57) 
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Integral i^ 3 -* can be written 


j( 3 ) — j( 3 ,l) _|_ 2 j( 3 ’ 2 ) _|_ /( 3 > 3 ) 

where 

r f r s 

/ (3 ’ 1} = / dsG(k,t,s) dr G(k,s,r)k 2 [ — G(k,t,r)]k 2 Q 
Jo Jo dk 

j(3,2) _ f ds G(k,t,s) [ drG(k,s,r)[k — G(k,t,r)][k 3 -^-] 

Jo Jo dk dk 

j(3,3) _ f ds G(k,t,s ) 2 ( dr G(k,s,r ) 2 k 4 —Q 
Jo Jo dk 

The term proportional to /( 3,1 ( can be written as 

J dk / (3,1) {(a (3) + a (6) )(WS + SW T ) + (d (3) + d {6) )(VU T S + SW)} D 
= -C( 7 V 7 ) - C( 8 V 8 ) + + 2a( 10 ) + 

where 

<7 (6) = s*- 0,0 -* {(a*- 3 -* + 2a( 6 ))[WS + ,SW T ] + (t/ (3) + 2d^)[VU T S + SVU]} D 

a (7) = s (bo){( a (3) + 2a (6) )[VUS + SVU T ] + (d (3) + 2d (6) )[VU T S + SVU}} D 

a (8) = s (o,D{( a (3) + 2a^)[VUS + SVU T ] + (d^ +2d^)[VU T S + SVU]} D 
<7 (9) = s (2,0) {(a (3) + 2a (6) )[WS + WR T ] + (d (3) +2d (6) )[W T S + ,SW]} D 

cr(io) = + 2a( 6 ))[WS + ,SW T ] + (t/ (3) + 2d^)[VU T S + SVU]} D 

<7 ( n) = s (0,2) {(a (3) + 2a (6) )[WS + WR T ] + (d (3) +2d (6) )[W T S + ,SW]} D 

and the quantities C^ 7 -* and C^ 8 -* are 

) (8) j 0 °°dk fcy 

Jq 00 dk kr 1 l 6 2 E 

Relaxation approximations for the time evolution of these terms are 

<7 (6) = T (3 ’ 0) - C'2 ,,O) 0<7 (6) 

& (7) = C (E°)Q a G) _ c£’ 0) ea (7) 

<7 (8) = T (3 ’ 3) - (58) 

b (9) = 2(7 (2 ’ o) 0cr (7) - c4 2 ,o) 0<7 (9) 

^(3°) = a (8) _ c a.i) 0(T (io) 

0(11) = T (3 ’ 2) - (7^ O ’ 2) 0cr (11) 
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To evaluate the /l 3,2 ) term, write 




where 


dk I (3 ’ 2) {(a (3) + a (6) )(VUS + SVU T ) + (d (3) + d (6) )(VU T S + SVU)} j 

= 3 a (r) +3a (S) _ a (12) _ a (13) 


a (12) = s (l, 0 ) « {kE) 


a (13) = 3 ( 0 , 1 ) « ( 


and the relaxation approximations are 

d (12) = -C'^ o,o), 0<7 (6) + 2C ( s 1 ’ oy Q(j (7) + - C ( k’ iY Q(J (12) 

where the modified transport coefficients are 


r,(o,o)'n f 0 °°dk(k V 'ffrE 

S f^dk 0jE S fWfaE 

r (..i)' n STdkWfkhE (I,,). 

s j™dkkri'\0 2 E R dk kr,'ie 2 (kEy 


where 


cr (13) = t (3 ’ 4) + 2T (3,2) - C , L O, 1 ), 0cj (13) 


(0;1) , / 0 °°dfc 2 r,kr,<\e 2 (kEy 

R / 0 °° dk kri'\e 2 (kEy 


To evaluate the contribution from /l 3 ’ 3 ), write 


& 4 <2" = 12k 2 Q - 8(k 3 Q)' + (& 4 <2)" 


Then 


dk / ( 3 ’ 3 ) {(a (3) + a (6) )(WS + ,SW T ) + (d (3) + d ( 6 ) )(W T S + SW)]^ 
= 12^) + a( 14 ) - 8a( 15 ) 


where 


all 4 ) = s^ 0,0 ) -j-(k 3 Q) 
dk 
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has the relaxation approximation 


and 


(7 


j( 14 ) = 2T (3,1) - c4 o,o), 0<t (14) 


a (15) =s (0,0)# k 4 Q 

dk 1 


(63) 


has the relaxation approximation 

d (15) = -c'g ) ’ o) "0r (3 ’ 1) + t (3,2) - c'2 ,,o),, 0<t (15) 

with 

W — fOO ;; L'fi W — roo 


f Q °° dk OlE 
To summarize, the second order model is 


f Q °° dk 0 1 (k 2 E)> 


(64) 

(65) 


- 3ct ( - 2 ' ) - + a (5) + 12(J (6) + (3 _ C ^)a^ + (3 - C™)*™ + 

+ 2cr ( - 10 ' ) + a^ 11 ) - a^ 12 ) - a^ 13 ) + a^ 14 ) - Sa^ 15 ) 

where the satisfy the system of relaxation equations Eqs. (55), (56), (58), (61), (63), 
(64) and the constants are defined by Eqs. (53), (54), (57), (60), (62), (65). 

Appendix II. Markovianized DIA 

Deriving useful conclusions from Leslie’s theory requires analytical expressions for the 
descriptors G^ 0 -* and Q of the isotropic background field. Suitable expressions will be 
obtained from a Markovianized DIA derived by evaluating the time dependence in the 
distant interaction limit. Let the real function iL(^), 0 < £ < oo satisfy 


H( 0) = 1, H(0< lfore>0, 


H(O d C < oo 


( 66 ) 


Then standard properties of delta functions imply 

— s)) ~ 8(t — s ) for A — > oo 


(67) 


Evaluate Eq. (10) in the distant interaction approximation in which A: — > 0 , p, q > oo. 
Then * 

7 (0) ( k ,M)= f dpdq B(k,p,q)G i0) (p,t,s)Q i0) (q,t,s) 

J k=p+q 

/ * f)n 

dp {k m ~7 — (k, p, p)G (0) (p, t, s)Q (0) (p, t, s ) 

UC[m 

~ B(K P, P)G (0) (P, t, s)k m p m p- ld Q—{p, t, s)} 


17 



where R(k,p,q) denotes the product of projection operators in Eq. (10) and f* denotes 
infrared regularization by restriction of the region of integration to (p, q) pairs satisfying 
p, q > k. 29 Assuming time stationary similarity forms G^°\p, t, s) = G^°\p r (t — s)) and 
Q(°\p,t, s ) = R(°\p r (t — s))Q^(p), the properties Eq. (66) of H may reasonably be 
postulated of the product G^ R^°\ Therefore, Eq. (67) implies that in this limit the 
damping is Markovian 

rj(°\]i,t,s) = 6(t — s)rj(°\ k) 

and Eq. (7) implies that the Green’s function is exponential, 

G^°\k, t, s) = exp [(s — f)?/°)(k)] for t > s (68) 

Likewise evaluating the force correlation Eqs. (11), (12) in the distant interaction limit 
implies that the forcing is white noise in time: 

</<“ > (k.()/ J (0| (k', s )>= <5(i-, s )<i(k + k')f’ ( " l (k) (69) 

Computing the correlation function from the relation 

Q (0) (k,t,s)(5(k + k') = f dn G w (k,t, ri ) [ dr 2 G (0) (k' ,t,r 2 )x 

Jo Jo 

</(°)(k,r 1 )/(°)(k',r 2 )> 

using Eqs. (68), (69) shows that the fluctuation dissipation relation 

Q (0 \k,t,s) = Q(°\k)[G w (k,t,s) + G (0 \k } s } t)] (70) 

expressing the time dependence of the correlation functions in terms of the response func- 
tion is also valid in this limit. 
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